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Abstract 

Predictive Higgs-fermion couplings can be obtained when a specific texture for 
the fermion mass matrices is included in the general two-Higgs doublet model. We 
derive the form of these couplings in the charged lepton sector using a Hermitian 
mass matrix Ansatz with four-texture zeros. The presence of unconstrained phases 
in the vertices (pihlj modifies the pattern of flavor-violating Higgs interactions. 
Bounds on the model parameters are obtained from present limits on rare lepton 
flavor violating processes, which could be extended further by the search for the 
decay r /i/i/x and /u — e conversion at future experiments. The signal from Higgs 
boson decays (pi TjjL could be searched at the large hadron collider (LHC), while 
e — jJL transitions could produce a detectable signal at a future e/i-collider, through 
the reaction ^~ ^ t~^t~ . 



PACS numbers: 12.60.Fr, 12. 15. Mm, 14.80. Cp 



1 Introduction. 



After many years of success of the Standard Model (SM), the Higgs mechanism is still 
the least tested sector, and the problem of electroweak symmetry breaking (EWSB) 
remmains almost as open as ever. However, the analysis of raditive corrections within 
the SM [1], points towards the existence of a light Higgs boson, which could be detected in 
the early stages of LHC [2] . On the other hand, the SM is often considered as an effective 
theory, valid up to an energy scale of 0{TeV), and eventually it will be replaced by a 
more fundamental theory, which will explain, among other things, the physics behind 
EWSB and perhaps even the origin of flavor. Several examples of candidate theories, 
which range from supersymmetry [3] to deconstruction [4], include a Higgs sector with 
two scalar doublets, which has a rich structure and predicts interesting phenomenology 
[5]. The general two-higgs doublet model (THDM) has a potential problem with flavor 
changing neutral currents (FCNC) mediated by the Higgs bosons, which arises when each 
quark type (u and d) is allowed to couple to both Higgs doublets, and FCNC could be 
induced at large rates that may jeopardize the model. The possible solutions to this 
problem of the THDM involve an assumption about the Yukawa structure of the model. 
To discuss them it is convenient to refer to the Yukawa lagrangian, which is written for 
the quarks fields as follows: 

£y = n^^L^iWi? + y2QL^2UR + Y,^QL<l>,dR + Y^Q^^^dR (1) 

where $12 = {(t>i2i 0i 2)"^ denote the Higgs doublets. The specific choices for the Yukawa 
matrices Y12 {q = u, d) define the versions of the THDM known as 1, 11 or 111, which 
involve the following mechanisms, that are aimed either to eliminate the otherwise un- 
bearable FCNC problem or at least to keep it under control, namely: 

1. DISCRETE SYMMETRIES. A discrete symmetry can be invoked to allow a given 
fermion type (u or d-quarks for instance) to couple to a single Higgs doublet, and 
in such case FCNC's are absent at tree-level. In particular, when a single Higgs 
field gives masses to both types of quarks (cither Y^ = Yf = or Y2 = Y2 = 0), 
the resulting model is referred as THDM-1. On the other hand, when each type of 
quark couples to a different Higgs doublet (either Y^ = Y2 = or Y2 = Yf = 0), 
the model is known as the THDM-II. This THDM-II pattern is highly motivated 
because it arises at tree-level in the minimal SUSY extension for the SM (MSSM) 
[5]. 

2. RADIATIVE SUPRESSION. When each fermion type couples to both Higgs dou- 
blets, FCNC could be kept under control if there exists a hierarchy between Y""''^ 
and y^"'^. Namely, a given set of Yukawa matrices is present at tree-level, but the 
other ones arise only as a radiative effect. This occurs for instance in the MSSM, 
where the type-II THDM structure is not protected by any symmetry, and is trans- 
formed into a type-HI THDM (see bellow), through the loop effects of sfermions 
and gauginos. Namely, the Yukawa couplings that are already present at tree-level 
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in the MSSM (Y-f, Y2) receive radiative corrections, while the terms {Y2, Y^) are 
induced at one-loop level. 

In particular, when the "seesaw" mechanism [6] is implemented in the MSSM to 
explain the observed neutrino masses [7, 8], lepton flavor violation (LFV) appears 
naturally in the right-handed neutrino sector, which is then communicated to the 
sleptons and from there to the charged leptons and Higgs sector. These corrections 
allow the neutral Higgs bosons to mediate LFV's, in particular it was found that the 
(Higgs- mediated) tau decay t ^ 3fi [9] as well as the (real) Higgs decay H Tfi 
[10], can enter into possible detection domain. Similar effects are known to arise in 
the quark sector, for instance B ^ /i/i can reach branching fractions at large tan/3, 
that can be probed at Run II of the Tevatron [11, 12]. 

3. FLAVOR SYMMETRIES. Suppression for FCNC can also be achived when a cer- 
tain form of the Yukawa matrices that reproduce the observed fermion masses and 
mixing angles is implemented in the model, which is then named as THDM-III. This 
could be done either by implementing the Frogart-Nielsen mechanism to generate 
the fermion mass hierarchies [13], or by studying a certain Ansatz for the fermion 
mass matrices [14]. The first proposal for the Higgs couplings along these lines was 
posed in [15, 16], it was based on the six-texture form of the mass matrices, namely: 

/ G \ 
Mi=\ CI Bi \. 
V Bt Ai I 

Then, by assuming that each Yukawa matrix Y12 has the same hierarchy, one finds: 
Ai ~ ms, Bi ~ ^777.2777.3 and Ci ~ ^Jm^rn^. Then, the Higgs-fermion couplings obey 
the following pattern: Hfifj ~ ^Jm^WTj / riiw ■, which is known as the Cheng-Sher 
Ansatz. This brings under control the FCNC problem, and it has been extensively 
studied in the literature to search for flavor-violating signals in the Higgs sector 
[17]. 

In this paper we are interested in studying the flavor symmetry option. However, the 
six-texture Ansatz seems disfavored by current data on the CKM mixing angles. More 

recently, mass matrices with four-texture Ansatz have been considered, and are found in 
better agreement with the observed data [18]. It is interesting then to investigate how 
the Cheng-Sher form of the Higgs-fermion couplings, gets modified when one replaces the 
six-texture matrices by the four-texture Ansatz. This paper is aimed precisely to study 
this question; we want to derive the form of the Higgs-fermion couplings and to discuss 
how and when the resulting predictions could be tested, both in rare tau decays and in 
the phenomenology of the Higgs bosons [10]. Unlike previous studies, we keep in our 
analysis the effect of the complex phases, which modify the FCNC Higgs couplings. 

The organization of the paper goes as follows: In section 2, we discuss the lagrangian 
for the THDM with the four-texture form for the mass matrices, and present the results 



3 



for the Higgs-fermion vertices in the charged lepton sector. Then, in section 3 we study 
the constraints impossed on the parameters of the model from low energy LFV processes. 
In section 4 we discuss predictions of the model for tau and Higgs decays, including the 
capabilities of future hadron and e/i-coUiders to probe this phenomena. Finally, section 
5 contains our conclusions. 



2 The THDM-III with four-texture mass matrices 



The Yukawa lagrangian of the THDM-III for the lepton sector is given by: 
After SSB the charged lepton mass matrix is given by. 



(2) 



Mi^-^{v,Yl + V2Y^), 



(3) 



We shall asuume that both Yukawa matrices Yj and Y-j have the four-texture form 
and Hermitic; following the conventions of [18], the lepton mass matrix is then written 
as: 



Ml 



Ci 
CI Bi Bi 
BI Ax 



when Bi one recovers the six-texture form. We also consider the hierarchy: 
I A; I I 5; I, I 5; I, I C; I, which is Supported by the observed fermion masses in the SM. 
Because of the hermicity condition, both Bi and Ai are real parameters, while the 

phases of Ci and Bi, ^b,c, can be removed from the mass matrix Mi by defining: Mi = 
P'^MP, where P = (iia(7[l, e**^, e**^*^+*'=')], and the mass matrix Mi includes only the 
real parts of Mi. The diagonalization of M is then obtained by an orthogonal matrix O, 
such that the diagonal mass matrix is: Mi = O^MiO. 

The lagrangian (2) can be expanded in terms of the mass-eigenstates for the neutral 
if^, A'^) and charged Higgs bosons {H"^). The interactions of the neutral Higgs bosons 
are given by. 



£1 = ^ 



g f m \ 



2 \mw / 



cos a 



Sij + 



+ 
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2 \mw/ 



cos j3 
sin a 



V2 sm{a - (5) f mw \ yi 
g cos /3 \ rrii ) ^'-^ 



cos (3 



+11 (^) J. 
2 \mw/ ' 



\/2 cos(a - (3) / mw \ yi 
g cos/3 V 

72 /mw\ y 



- tan pSij + ( ^ ) r^. 

g cos P \ nii J ■' 



rrii ) 



(4) 
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The first term, proportional to 6ij corresponds to the modification of the THDM-II over 
the SM result, while the term proportional to Yj denotes the new contribution from 
THDM-III. Thus, the fermion-Higgs couphngs respect CP-invariance, despite the fact 
that the Yukawa matrices include complex phases; this follows because of the Hermiticity 
conditions impossed on both Y[ and 1^2- 

The corrections to the Icpton fiavor conserving (LFC) and fiavor-violating (LFV) 
couplings, depend on the rotated matrix: Vj = O^PYjP^O. We shall evaluate Vj, by 
assuming that Y2 has a four-texture form, namely: 



Cj 
C2 B2 
B* 



^2 I > I ^2 



Bo 



Co 



(5) 



The matrix that diagonalizes the real matrix Mi with the four-texture form, is given 



by: 



O 



A2A3(A-Ai) 
A(A2-Ai)(A3-Ai) 



Ai(Ai-A) 
(A2-Ai)(A3-Ai) 

Ai(^-A2)(^-A3) 
A(A2-Ai)(A3-Ai) 



AiA3(A2-A) 
A(A2-Ai)(A3-A2) 

A2(A-A2) ~ 
(A2-Ai)(A3-A2) 

\2{A-\l){\3-A) 

A(A2-Ai)(A3-A2) 



AiA2(A-A3) \ 
A(A3-Ai)(A3-A2) 

A3(A3-A)~ 
(A3-Ai)(A3-A2) 

A3(A-Ai)(A-A2) 
A(A3-Ai)(A3-A2) J 



where Trip 



mi 



Ai \,m. 



1712 



A. 



A3 |,?7 = Aa/ma 



Then the rotated form Y2 has the general form. 



= O^PY^P^O 

^J211 ^J212 -'_213 

^221 ^222 ^223 | • (6) 
^ 231 ^ 232 ^ 233 

However, the full expressions for the resulting elements have a complicated form, as 
it can be appreciated, for instance, by looking at the element (^^2)22, which is displayed 
here: 



777,3 - A2 V ^"^2 "73 - A2 




(7) 



+A,^^ - \Bie»' + fl.e-*-]. '^-''"7-^' (8) 

7773 - A2 V ?773 - A2 

where we have taken the limits: 1^41 , 771,-, rn^ ^ rrie- The free-parameters are: -B2, -B2, ^2, ^• 
To derive a better suited approximation, we shall consider the elements of the Yukawa 
matrix Y2 as having the same hierarchy as the full mass matrix, namely: 
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C2 = C2\l J (9) 



B2 = b2^JiA-X2)im3-A) (10) 

B2 = 62(^3 -A + A2) (11) 

A2 = 02^. (12) 

Then, in order to keep the same hierarchy for the elements of the mass matrix, we 

find that A must fall within the interval (ma — 7712) < A < m^. Thus, we propose the 
following relation for A: 

A^m^{l- (3z), (13) 

where z — 777,2/777,3 -C 1 and < /3 < 1. 

Then, we introduce the matrix x as follows: 



Xij 



ij V 

n ■ 

"Xije"- (14) 



which differs from the usual Chcng-Shcr Ansatz not only because of the appearence of 
the complex phases, but also in the form of the real parts Xij = \Xij\- 

Expanding in powers of z, one finds that the elements of the matrix % have the 
following general expressions: 

Xii = [fe2-(c;e^*^+C2e-^*^)]77+[a2 + 62-(&;e^*"+&2e-^*")]/? 

X12 = (c2e-**^-62)-r7[a2 + 62-(6;e**^+62e-^*^)]/?] 

Xi3 = ia2-b2e-'''^)ri^ 

X22 = 6277+[a2 + 62-(fe;e'*" + 62e-**")]/3 

X23 = (fe2e-^*--a2)V^ 

X33 = 0,2 (15) 

It is also relevant to point out the following: 

• When the phases $b and vanish, [3=1 and one takes the 6-texture limit 
{B2 0, i.e. b ^ ^ 7] = -1), Eq. (14) reduces to 



(Yj)^^ = (2c2 + a2-2b2)mi/v 
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= (C2 + 02 - 262) V"^l"^2/ 




= (62 - 02) y/miTTls/v 




= {a2-2b2)m2/v 




= (62 - 02) \/m2m3/v 




= a2mz/v 



which correspond to the Unsafe of Cheng-Sher (See Eq. (32) in Ref. [15]). 

• On the other hand, when the phases $b and $c vanish, /3 — 777.2/7713 and rj = 1, 
Eq. (14) reduces to 



(^2)11 = (62-2c2)mi/^; 

(■^2)12 ^ (c2 - &2) y/mim2/v 

(^2)^3 = (02 - &2) Vrnlrni/v 

(^2)22 = ^2-2/- 

(^2)23 ^ (62-02)^2/^ 

(y2033 = «2m3A (17) 
in this case one reproduces the results given in Ref. [20] (See Eq. (24)). 



While the diagonal elements xu are real, we notice (Eqs. 15) the appearance of 
the phases in the off-diagonal elements, which are essentially unconstrained by present 
low-energy phenomena. As we will see next, these phases modify the pattern of flavor 
violation in the Higgs sector. For instance, while the Cheng-Sher Ansatz predicts that 
the LEV couplings (i^')i3 and {¥2)23 vanish when 02 = &2) in our case this is no longer 
valid for cos$b 7^ 1. Furthermore the LEV couplings satisfy several relations, such as: 
1X23 1 = IXisL which simplifies the parameter freedom. 

Finally, in order to perform our phenomenological study we find convenient to rewrite 
the lagrangian given in Eq. (4) in terms of the x^'s as follows: 



/ mj \ cos a ^ ^ sin(Q; - P) f y^mirrij \ ^ 
\ mw ) cos 13 V2 cos V "^vk / 



-I — — 1 



mw) cos/5 1/2 cos/5 

- ( J tan /5 



1 / y/miiTij 
"\/2 cos /5 \ mw 



mi \ sin a cos(a - (3) ( ^m^m^ \ ^ 

I ~ ^ii H ~ I I Xij 



Xij 



(18) 
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where, unlike the Cheng-Sher Ansatz, Xij {i ^ j) are complex. 



3 Bounds on the LFV Higgs parameters 

Constrains on the LFV-Higgs interaction will be obtained by studying LFV transitions, 
which include the 3-body modes {k — > Ijlklk), radiative decays — > Ij + /i — e 
conversion in nuclei, as well as the (LFC) muon anomalous magnetic moment. 

3.1 LFV three-body decays. To evaluate the LFV leptonic couplings, we calculate the 
decays li — > IjlJk: including the contribution from the three Higgs bosons and 
A^).We obtain the following expression for the branching ratio: 

„ ,, ,,Tv 5(^,-fc + 2 Ti mnrriim^ { cos^ ia — 13) siv? a sin^fo; — /3) cos^ a 
Brik Llkk) = — TT-^ — — 7 — - { \ 

2 cos(q; — (3) sin(Q; — (5) cosa sin a sin^ (3 1 xfj 

mlp mjjo m\o J 2 cos^ /? 

where Tj denotes the life time of the lepton 1^ and we have assumed Xkk <^ 1; this result 
agrees with Ref. [20]. 

In particular, for the decay r~ — > ii~ fi'^ ii~ we obtain the following expression for the 
branching ratio: 

„ , _ _ . h Tr m%m% fcos^fo; — /3) sin^ q; siv?' ia — (3) cos^ a 



2/3^ ,.2 



cos(q; — f3) sin(a — /3) cos a sin a sin /3 1 X23 /^,^^ 
mlomjjo m\o J cos*/? 

here r^- corresponds to the hfe time of the r lepton (we have also assumed X22 <S 1)- 

Using the experimental result Br{T~ — > /i^ /i'^ iJ,~) < 1.9 x 10~^, we get an upper 
bound on X23 {{X23)u~b^^) as a function of a and tan/?. In Fig. 1 we show the value of 
this bound as a function of tan/5 for a = /3 — 7r/4, a = /3 — 7r/3 and a = (3 — 7r/2, taking 
m/jO = 115 GeV and uiho = m^o = 300 GeV. 

Taking X23 ~ 1, tan/? f^i 30 and 7r/4 < /? — a < 7r/2, in Eq. (20) one finds tipically 
that Br{T — > 3//) ~ 10~^, which puts it into the regime that is experimentally accessible 
at r-factorics over the next few years. At LHC and SuperKEKB, limits in the range of 
10~^ should be achievable [21], allowing a deeper probe into the parameter space. 

3.2 Radiative decays. The branching ratio of /i'^ — > e'^j at one loop level is given by [22] 



D / + + . aemTumimimi 2 2 



cos^(q; — 13) 



m 



1 ^4 
In^ 



3 

+ - 



8 



+2 



cos^(o; — /?) sin^(Q; — (3) 



2 2 



+ 



sm 



{a -(3) 



m 



HO 



In 



In 



3 

+ 2 



+ 



3 

+ 2 



In 



AO 



In 



3 

+ 2 



3 

+ - 



m 



AO 



(21) 



(a2 — 626 ) I y/p. We will make use of the current 
e+7) 



Prom Eqs. (15) we have X23 = X13 

experimental upper bound Br{jjL^ — > e'^7) < 1.2 x lO"-*^^ [23] to constraint X23(xi3) as 
a function of a and tan/?. Assuming m/jO = 115 GeK and iriijo = m^o = 300 GeV, we 
depict in Fig. 2 the value of the upper bound on X23 {{X23)u^'^"') as a function of tan/?, 
again for a = /? — 7r/4, a = /? — 7r/3 and a = P — 7t/2. A new experiment at PSI will 



measure the process — > e+7 with a sensitivity of 1 event for Br{iX^ e+7) = 10" 
[24], which would improve the upper bound on X22, by a factor ~ 10~^/^ 0.18 . 
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3.3 n ~ e conversion. The formulas of the conversion branching ratios for the Icpton 
flavor violating muon electron process in nuclei at large tan jS, in the aluminum and lead 
targets, are approximately given by 



Br{fj,-Al e-Al) 



. mi m, ml tan® 3 cos^ 3 „ 
1.8 X 10-^ '"4 

2 m]jo IjJcapt 



X12 



and 



Br{ii-Pb e-Ph) ~ 2.5 x 10 



3 mi 777.2 "^p ^^^^ COS^ (5 2 

Xl2' 



2 ti^ m|^o iOcapt 



(22) 



(23) 



respectively, where uJcapt is the rate for muon capture in the nuclei [25]. The values are 
<^copt — 0.7054 X 10^ and oJcapt — 13.45 x 10® in the aluminum and the lead nuclei, 
respectively [26]. There are several planned experiments which are aiming at improving 
the bounds of the branching fractions for relevant processes by three or four orders of 
magnitude [27, 28, 29] . In particular, the MECO experiment will search for the coherent 
conversion of muons to electrons in the field of a nucleus with a sensitivity of 1 event for 
5 X 10^® muon captures, i. e. Briii'M e'M) < 2 x 10"^^ [30, 31]. Taking mj^o = 300 
GeV, we plot in Pig. 3 the value of the upper bound on xu {{Xi2)u 



as a function 



of tan/9 for Al and Pb, for the current experimental measurement Br{fi M 
6.1 x 10~^^ [31]. In Pig. 4, we show the same as in Pig. 3 but for Br{ii~J\f 
2 X 10-^^. 



e-M) < 
e-M) < 



3.4 Muon anomalous magnetic moment. Taking the average value of the measurements 
of the muon {g — 2) from [33] and the recent analysis by different groups [34, 35] one can 
conlude that 

Aa^ = a^^P - a^^ 300 ± 100 X 10-^^ (24) 

The contribution to the muon — 2 of the one loop level fiavor changing diagram is 
given as follows. 
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Aa 



± 



1 m\m\ cos'(('\ — /3) 



( 



In 




(25) 



where the sign +(-) is for scalar (jp — (pseudoscalar, 0° = A°) exchanges [32, 36, 37, 
38]. From Eq. (24) is clear that we need to increase the theoretical value of a^. Hence, 
we will consider the contribution of to the muon {g — 2) assuming X23 = 1, namely 
Aa^ (X23 = !)• We take m^o = 115 GeV and present in Fig. 5 the result for such 
contribution as a function of tan/3 for a = /3 — 7r/4 and a = (3 — tt/S. We observe that 
Aa|j°(X23 = 1) < 240 X 10~^^. On the other hand, the contribution to the anomalous 
magnetic moment from two-loop double scalar-exchanging diagrams is comparable with 
the one from the corresponding flavor- changing one loop diagrams [32]. It was already 
shown in Ref. [32] that the two-loop double scalar (pseudoscalar) exchanging diagrams 
give negative (positive) contributions, which have opposite signs as the one from one loop 
scalar (pseudoscalar) exchanging diagram. Hence, we can conclude that it would be very 
hard to constrain X23 from the muon g — 2 measurements, or to explain such deviation 
from the pure Higgs sector in case the signal is confirmed. 

Thus, we conclude from this section that the bounds on the LFV parameters are given 
as follows. 

• X12 < 5 X 10"^, from /i~ — e~ conversion experiments. 

• Xi3 — X23 < 6 X 10~^, from the radiative decay e+7 measurements. 

However, one can still say that at the present time the couplings Xi/^ highly con- 

strained, thus they could induce interesting direct LFV Higgs signals at future colliders. 

4 Probing the LFV Higgs couplings at future collid- 
ers 

In order to probe the LFV Higgs vertices we shall consider both the search for the LFV 
Higgs decays at future hadron colliders (LHC mainly), as well as the production of Higgs 
bosons in the collisions of electrons and muons, which was proposed some time ago [39], 
namely we shall evaluate the reaction e/x — > /i° ^ rr. 

4.1 Search for LFV Higgs decays at Hadron colliders. We shall concentrate here on the 
LFV Higgs decays 0j t/j,, which has a very small branching ratio within the context of 
the SM with light neutrinos (< 10~^ — 10~^ ), so that this channel becomes an excellent 
window for probing new physics [10, 40, 41]. The decay width for the procces 0i Tji 
(adding both final states t'^^~ and t~ ^'^ ) can be written in terms of the decay width 
r{Hi ^ tt), as follows: 



r(0, ^ r//) = {R^S 



nH, ^ tt) 



(26) 
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where 



^'^ 94>TT COS a y ruj 



Therefore, the Higgs branching ration can be approximated as: Br{(j)i t/j.) = (Rf^i) x 



Sr(0i tt). We calculated the branching fraction for h r//, and find that it reaches 
values of order 10~^ in the THDM-III; for comparison, wc notice that in the MSSM case, 
even for large values of tan/?, one only gets Br{h — > Tji) ~ 10^'^. 

These values of the branching ratio enter into the domain of detectability at hadron 
colliders (LHC) , provided that the cross-section for Higgs production were of order of the 
SM one. Large values of tan (3 are also associated with large b-quark Yukawa coupling, 
which in turn can produce and enhancement on the Higgs production cross-sections at 
hadron colliders, even for the heavier states and cither by gluon fusion or in the 
associated production of the Higgs with b-quark pairs; some values are shown in table 1; 
these were obtained using HIGLU [42] . Thus, even the heavy Higgs bosons of the model 
could be detected through this LFV mode. 



mH,A [GeV] 


<a [Pb] 


< [pb] 


< [pb] (^ CX-) 


150 


126.4 (492.6) 


129.1 (525.) 


200 (800) 


200 


29.5 (114.3) 


29.1 (120.) 


100 (400) 


300 


3.6 (13.5) 


3.15 (13.6) 


20 (80) 


350 


1.6 (5.9) 


1.2 (5.6) 


12 (48) 


400 


0.75 (2.75) 


0.73 (2.8) 


8 (32) 



Table 1. Cross-section for Higgs production at LHC, through gluon fusion {cr^^^) and in 
association with bb quarks, (cr^'"^), for tan/3 = 30 (60). 

For instance, for mH,A = 150 GeV and tan/3 = 30(60) the cross-section through gluon 
fusion at LHC is about 126.4 (492.6) pb [42], then with Br{H rfj,) ~ 10-^(10-^) and 
an integrated luminosity of 10^ pb~^, LHC can produce about 10'^(10'') LFV Higgs events. 
In Ref. [43] it was proposed a series of cuts to reconstruct the hadronic and electronic tau 
decays from h ^ Tfi and separate the signal from the backgrounds, which arc dominated 
by Drell-Yan tau pair and WW pair production. According to these studies [43], even 
SM-like cross sections and ~ 150 GeV, one coud detect at LHC the LFV Higgs decays 
with a branching ratio of order 8 x 10~^, which means that our signal is clearly detectable. 

4.2 Tests of LFV Higgs couplings at e/i- colliders. Another option to search for LFV Higgs 
couplings, but now involving the electron-muon-Higgs couplings, would be to search for 
the reaction: e~(pa) + IJ''^(pb) — > — > T"~{pc) +T'^{Pd)- Assuming X33 <^ 1, the result 
for the cross section is given by: 

, _ , S 17111712 171^ 2 

ale Li ^ T T ) — — -. rTTXi9 

^ ^ ^ 32 7rv4 cos4/3-^^2 
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{|D/io(s) p cos^(q; — /3) sin^ a 

—2Re{Di^o{s) D^o(s)} cos(q; — P) sm{a — [3) sin a cos a 

+ 1 Dho (s) 1 2 sin^ {a- (5) cos^ a + | D^o (s) | ^ sin^ /?} , (28) 

where D^o{s) denotes the Breit-Wigner form of the 0° propagator 

D^o{s) = {s- mjo + im^o r|2)"^ (29) 

and s^{pa+ PhY = (Pc + Pd)^- 

The non-observation of at least an event in a year would imply that 

(T(e~/x^ — > T^r"^) X luminosity x lyear < 1, (30) 

which would allow us to put an upper bound on Xi2^ namely {xi2Tub''^'^{s) as a function 
of a and tan/3. In order to obtain numerical results, we take T^l^ = 0.004 GeV for rriho — 
115 GeV; Fg^ = 0.14 GeV for mno = 300 GeV; and rf, = 0.045 GeV for m^o = 300 GeV 
[44] and a luminosity £ = 2 x 10^^ cm~^ [39]. We present our numerical results for 
(Xi2)rr^^(s = iTT'lo) and {xi2)Tb!^'^ i^'^ = '^ho = "^^o) in Fig. 6 and Fig. 7, respectively. 

We can also estimate the number of events N'^^~*'^'^ (s) by taking for Xi2 the value for 
the current upper bound on Xi2 obtained from the measurements in /x^ — conversion 
experiments, namely Xi2 — {Xi2)'^^^ as a function of tan/3 for Al, for the current 
experimental measurement Br{ijL~J\f ^ e~J\f) < 6.1 x 10~^^ [31]. Hence, we get 

N"'^^'^'' (s) = aie'ii" T^T^) X luminosity x lyear, (31) 

as a function of a and tan/3. In order to obtain numerical results, we take T^°^ — 
O.OOAGeV for m^^o = lib GeV; Fg" = 0.14 GeV for m^o = 300 GeV; and F^^^ = 
0.045 Gel/ for uiao = 300 GeV [44] and a luminosity £ = 2 x 10^2 cm^^ g-i [39]. We 
present our numerical results for N^^^~''^'^{s = m^o) and N'^'^"*'^'^ [s = mjjo = "^io) in Fig. 
8 and Fig. 9, respectively. We obtain around 100 events per year, which is very likely 
detectable. 



5 Conclusions 

We have studied in this paper the lepton-Higgs couplings that arise in the THDM-Ill, 
using a Hermitic four-texture form for the leptonic Yukawa matrix. Because of this, 
although the fermion-Higgs couplings are complex, the CP- properties of (even) 
and A'^ (odd) remmain valid. 

We have derived bounds on the LFV parameters of the model, using current experi- 
mental bounds on LFV transitions. Our resulting bounds can be summarized as follows. 

• X12 < 5 X 10~^, from ii~ — e" conversion experiments. 

• Xi3 — X23 < 6 X 10"^, from the radiative decay — > measurements. 
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However, one can say that the present bounds on the couphngs Xy's still allow the 
possibility to study interesting direct LFV Higgs signals at future colliders. 

In particular, the LFV couphngs of the neutral Higgs bosons, can lead to new discovery 
signatures of the Higgs boson itself. For instance, the branching fraction for H/A — > r/x 
can be as large as 10~^, while Br{h r/i) is also about 10^^. These LFV Higgs modes 
complement the modes fifi, r 3/i, t ^ fi'-f and fi 67, as probes of flavor 

violation in the THDM-III, which could provide key insights into the form of the Yukawa 
mass matrix. 

On the other hand, one can also relate our results with the SUSY-induced THDM-HI, 
by considering the effective Lagrangian for the couplings of the charged leptons to the 
neutral Higgs fields, namely: 

-£ = LlYM^, + LlYi (ei 1 + es^Jy.) Ir4* + h.c. (32) 

In this language, LFV results from our inability to simultaneously diagonalize the term 
Yi and the non-holomorphic loop corrections, e2YiYjYiy. Thus, since the charged lepton 
masses cannot be diagonalized in the same basis as their Higgs couplings, this will allow 
neutral Higgs bosons to mediate LFV processes with rates proportional to e|. In terms 
of our previous notation we have: Y2 = e2YiYjY,y. An study of the values for €2 resulting 
from general soft breaking terms in the MSSM is underway. 
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Figure Captions 

Fig. 1: The upper bound {X23)u^^'^ as a function of tan /3 for a = /3 — 7r/4, a — P — 7r/3, 
a = /3 - 7r/2, with Br{T~ /i" /i'^ /i") < 1.9 x 10"^, taking niho = 115 GeV and 
mjjo = m^o = 300 GeV. 

Fig. 2: The upper bound {X23)u^^'^ as a function of tan /3 for a = /3 — 7r/4, a — P — tt/S, 
q; = /3 - 7r/2, with Br{iJ& ^ e+7) < 1.2 x 10~^^, taking m/jo = 115 GeF and niHO = 
m^o = 300 GeV. 

Fig. 3: The upper bound (xi2)l!^r''^ as a function of tan (5 for A/, P6, with Br{fjL-J\f 
e N) < 6.1 X 10"^^ and assuming uiho ^ 300 GeV. 

Fig. 4: The same as in Fig. 3, but taking Br{iJi~M e~M) < 2 x 10"^'''. 

Fig. 5: Aa|j° as a function of tan/3 lov a — (3 — tt/A, a — (3 — — 7r/3, with X22. — 1 and 
assuming m/jo = 115 GeV . 

Fig. 6: The upper bound {xi2)Zb^^^ for « = "^^0 = (115G'ey)2, with V^l^ = 0.004^61/, 
as afunction of tan /3 for a = — 7r/4, a = /3 — 7r/3, when (T(e~//^ — > r^r"^) xluminosity x 
1 year < 1, taking £ = 2 x 10^^ cm"^ s"-^ 

Fig. 7: The upper bound {xi2)Zb^^^ for s - m^^o = m^o = (300^6^)2, with = 
0.14 Gey and = 0.045 GeF, as afunction of tan /3 for a = /3-7r/4, o; = /3-7r/3, a = 
/3 — 7r/2, when (7(e~//+ ^ t~t+) x luminosity x 1 year < 1, taking C — 2x 10^^ cm~^ s"-^ 

Fig. 8: Number of events N^t"^^^ for s = m^o = {llbGeVf, taking X12 = (Xi2)(^l'r''^' 
for Br{n-Al e'Al) < 6.1 x 10-^^ with r'l°^ = 0.004 GeF, as a function of tan/3 for 
a = /3 - 7r/4, a = /3 - 7r/3, taking £ = 2 x 10^^ ^^-2 ^-i 

Fig. 9: Number of events N^^^^'^'^ for s = mjjo = rn^o = (300 GeV^, taking X12 — 
(Xi2)(^l'r'^' for Br{^-Al AV) < 6.1 x 10~^^ with TfJt = 0.14 GeV and t£ = 

0.045 GeV, as a function of tan/3 ior a — P — 7r/4, a — P — 7r/3, a — P — 7r/2, taking 
£ = 2 X 10=^2 g-i 
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